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Plasmon  polaritons  in  photonic  superlattices  containing  a left-handed  material 
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We  analyze  one-dimensional  photonic  superlattices  which  alternate  layers  of  air  and  a left-handed 
material.  We  assume  Drude-type  dispersive  responses  for  the  dielectric  permittivity  and  magnetic 
permeability  of  the  left-handed  material.  Maxwell’s  equations  and  the  transfer-matrix  technique 
are  used  to  derive  the  dispersion  relation  for  the  propagation  of  obliquely  incident  optical  fields. 

The  photonic  dispersion  indicates  that  the  growth-direction  component  of  the  electric  (or  magnetic) 
field  leads  to  the  propagation  of  electric  (or  magnetic)  plasmon  polaritons,  for  either  TE  or  TM 
configurations.  Furthermore,  we  show  that  if  the  plasma  frequency  is  chosen  within  the  photonic 
(n(u)))  = 0 zeroth-order  bandgap,  the  coupling  of  light  with  plasmons  weakens  considerably.  As  light 
propagation  is  forbidden  in  that  particular  frequency  region,  the  plasmon-polariton  mode  reduces 
to  a pure  plasmon  mode. 
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Over  the  years,  artificial  complex  materials  have  been 
increasingly  used  to  shape  and  manipulate  light—————. 
The  microstructuring  of  high  quality  optical  materials 
yields  remarkable  flexibility  in  the  fabrication  of  nanos- 
tructures, and  allows  for  the  tailoring  of  electromagnetic 
dispersions  and  mode  structures  to  suit  almost  any  need. 
Metamaterials&Zi8!9!10!11:12]13,  also  known  as  left-handed 
materials  (LHMs),  are  a remarkable  example  of  such 
nanostructuring.  Light  propagation  through  metamate- 
rials is  characterized  by  a phase  velocity  opposite  to  the 
Poyinting  vector,  which  corresponds  to  negative  disper- 
sive electric  and  magnetic  responses. 

The  advent  of  metamaterials  has  opened  a new  era  for 
optical  devices,  and  has  also  given  considerable  thrust 
to  the  recent  area  of  plasmonics,  the  study  of  plasmon 
polaritons.  In  metal-dielectric  interfaces,  for  example, 
surface-plasmon  polaritons  are  coupled  modes  that  re- 
sult from  resonant  interactions  between  electromagnetic 
waves  and  mobile  electrons  at  the  surface  of  a metal  or 
semiconductor.  Such  resonant  surface-plasmon  polari- 
tons may  have  much  shorter  wavelengths  than  that  of 
the  radiation,  which  enables  them  to  propagate  along 
nanoscale  systemsi^i&i,  opening  up  a wide  range  of  pos- 
sibilities for  the  construction  of  new  optical  devices.  It  is 
well  known,  for  example,  that  one  of  the  most  important 
features  of  surface-plasmon  polaritons  is  to  confine  light 
to  very  small  dimensions,  yielding  the  merging  of  pho- 
tonics and  electronics  at  the  nanoscale.  Furthermore, 
considering  the  dispersive  character  of  the  LHM’s,  to- 
gether with  the  enhanced  optical  magnetism  exhibited 
by  them,  one  might  conjecture  the  existence  of  remark- 
able new  phenomena  such  as  the  excitation  of  plasmon 
polaritons  of  a magnetic  nature,  that  is,  magnetic  den- 
sity waves  resulting  from  resonant  interactions  between 
the  optical  field  and  current  densities  at  the  metamate- 


rial. Ultimately,  those  developments  could  lead  to  in- 
creases in  the  resolution  of  microscopes,  in  the  efficiency 
of  LEDs,  and  in  the  sensitivity  of  chemical  and  biological 
devices^&i^. 

One-dimensional  (ID)  superlattices  which  al- 
ternate layers  of  positive  and  negative  mate- 
rials have  already  exhibited  many  interesting 
properties14i15i16i17i18i19i20’21i22i23  that  are  absent  in 
superlattices  composed  solely  of  positive  materials. 
In  particular,  the  existence  of  a non-Bragg  photonic 
bandgap,  also  known  as  a zeroth  order  gap,  has  been 
suggested^,  detected13'14,  and  characterized15!16.  In 
order  to  investigate  the  possibility  of  excitation  of 
electric/magnetic  plasmon  polaritons,  in  the  present 
work  we  study  the  oblique  incidence  of  light  on  a model 
ID  superlattice  composed  of  layers  A of  air,  and  layers 
B of  a doubly  negative  material.  Layers  A (width  a)  and 
B (width  b)  are  distributed  periodically  so  that  d = a + b 
is  the  period  of  the  superlattice  nanostructure  (cf.  Fig. 

HD. 

In  the  B layers,  the  electric  and  magnetic  responses 
are  dispersive  and  may  assume  negative  values.  If  one 
neglects  losses,  they  may  be  described  by^i— iM -id 

eB{u>)=e0-'4  ! = (!) 

where  £_b(w)  and  hb{^)  are  the  dielectric  permittivity 
and  magnetic  permeability  in  slab  B,  respectively,  one 
may  chooseiiii  £o  = 1.21  and  no  = 1.0,  and  the  elec- 
tric/magnetic plasmon  modes  are  at  v = ve  = anc^ 

v = vm  = which  correspond  to  the  solutions  of 

£b(w)  = 0 and  hb{w)  = 0,  respectively. 

We  note  that  dispersions  such  as  those  in  JT])  hold 
in  periodically  LC  loaded  transmission  linesii.  These 
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negative-index  systems  were  shown  to  exhibit  good  mi- 
crowave properties,  with  low  loss  and  broad  bandwidthAA. 

We  shall  be  interested  in  studying  the  properties  of 
both  the  transverse-electric  (TE:  electric  field  parallel  to 
the  interface  plane,  see  Fig.  1)  and  transverse-magnetic 
(TM:  magnetic  field  parallel  to  the  interface)  polariza- 
tions of  a monochromatic  electromagnetic  field  of  fre- 
quency to  propagating  through  a ID  periodic  system. 

In  the  case  of  a TE  field,  one  may  choose 

E(r,  t)  = E(z)  exp  [i  (qx  — ust)]  ey  , (2) 

whereas,  in  the  case  of  TM  polarization,  the  magnetic 
field  may  be  considered  as 


H(r  ,t)  = H(z)  exp  [i  (qx  — u>t)]  ey,  (3) 


where  we  have  assumed  that  the  superlattice  was  grown 
along  the  z axis,  q is  the  wavevector  component  along  the 
x direction,  and  ey  is  the  cartesian  unitary  vector  along 
the  y direction. 

Maxwell’s  equations  lead  to  the  following  differential 
equations  for  the  electric  and  magnetic  amplitudes 
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where  n(z)  = -J z(z)\J y,(z)  is  the  ^-position  dependent 
refraction  index. 

In  the  sequel,  n A and  nB  are  the  refraction  indices, 
while  and  Mb  are  the  magnetic  permeabilities  in  lay- 
ers A and  B,  respectively.  Note  that,  because  of  CD.  the 
ub  refraction  index  in  layer  B is  a function  of  the  fre- 
quency, and  may  be  a real  positive,  real  negative  or  a 
pure  imaginary  number.  Also,  note  that  q,  the  incident 
wavevector  component  along  the  x direction,  may  be  ob- 
tained as  a function  of  the  angle  of  incidence  9 = 6a  by 
q = ~nA  sin#. 

Equations  © and  © may  be  solved  by  means  of  the 
transfer-matrix  technique  (see,  for  example,  the  work  by 
Cavalcanti  et  al — and  references  therein).  For  nB  a real 
number,  and  n2B  — ei2a  sin2  9 > 0,  the  procedure  yields  the 
TE  polarization  dispersion  relation  from  the  solution  of 
the  transcendental  equation 


cos (kd)  = cos(Q Aa)  cos(Q Bb) 

~ (|0  + s'm(QAa)sin(QBb).  (6) 

In  the  formula,  k is  the  Bloch  wavevector  along  the  z di- 
rection which  is  the  axis  of  the  photonic  crystal;  fields  in 
consecutive  unit  cells  are  related  by  the  Bloch  condition, 
i.e.,  by  the  phase  factor  elkd.  Qa  and  Qb  are  defined  as 
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FIG.  1:  (Color  online)  Pictorial  view  of  the  ID  multilayer 
photonic  superlattice  with  layers  A and  B in  periodic  arrange- 
ment, and  the  electric  and  magnetic  fields  for  the  TE-like  and 
TM-like  electromagnetic  waves  schematically  shown.  Note 
that,  for  normal  incidence,  the  two  polarizations  are  equiva- 
lent. 
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FIG.  2:  (Color  online)  TE  and  TM  dispersion  relations  v = 

u(k)  in  photonic  periodic  superlattices  (is  = ■£-).  Calculations 

were  performed  for  air  as  slab  A (ca  = 1,  yA  = 1),  a = b = 12 

mm  (also  for  a = 12  mm  and  b = 6 mm),  and  uje/2n  = 

uJm / 27t  = 3 GHz  for  the  Drude  model  [cf.  Eq.  (©]  in  slab 

B.  Dotted  lines  indicate  the  pure  electric/magnetic  plasmon 

modes  at  v = ve  = „ — and  v = Vm  = gm_.  Solid  and 
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dashed  lines  correspond  to  9 = 7t/12  and  9 = 0 incidence 
angles,  respectively. 


and 

Qb  = ~.\ln2B  ~n2Asm29  = ^-\nB\\cos9B\,  (8) 

where,  FA  = (Qa/ha),  Fb  = (Qb/vb)  and,  in  the 
last  equality,  we  have  made  use  of  Snell’s  law  ha  sin  9 = 
nB  sin  9b- 

For  nB  a real  number,  and  n2B  — ri2A  sin2  9 < 0,  the 
procedure  yields  the  TE  polarization  dispersion  relation 
from  the  solution  of  the  transcendental  equation 

cos  (kd)  = cos(Q  Aa)  cosh(Q  Bb) 

~\  “ 07)  sin(<2Aa)sinh(Qs&).(9) 

where  Qa  is  still  given  by  0,  but 


Qa  = —ua | cos  6\  = —nA \ cos  0A\,  (7) 
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FIG.  3:  (Color  online)  TE  dispersion  relations  v = v(k ) in 
photonic  periodic  superlattices  {v  = as  functions  of  vari- 
ous incidence  angles.  Calculations  were  performed  using  the 
same  parameters  as  in  Fig.  2.  Dotted  lines  indicate  the 
pure  magnetic  plasmon  modes  at  v = vm  = 2^^L—,  dashed 
lines  correspond  to  0 = 0,  and  solid  lines  correspond  to  (a) 
0 = 7t/12,  (b)  0 = 7t/6  and  (c)  0 = 7r/3  incidence  angles, 
respectively. 


FIG.  4:  (Color  online)  TE  and  TM  dispersions  in  a photonic 
periodic  superlattice  with  a = b = 12  mm.  Numerical  calcu- 
lations were  carried  out  for  air  as  slab  A (ea  = 1 and  (ia  = 1). 
For  slab  B,  eb  and  /rs  were  chosen  by  taking  plasma  frequen- 
cies in  the  Drude  model  as  u}e/2n  = 1.5  GHz  and  ajm/27r  = 2.5 
GHz,  respectively,  in  (a)  and  (b),  whereas  in  (c)  and  (d)  we 
have  taken  we/27r  = 2.5  GHz  and  uim/2n  = 1.5  GHz.  Note 
that,  in  panel  (b),  the  electric  plasmon  frequency  i/e  = 2J^_ 
is  in  the  (n(u>))  = 0 zeroth-order  bandgap  region,  leading 
to  a basically  flat  0 = 7r/12  plasmon-polariton  band  with  a 
frequency  essentially  equal  to  ve  for  the  TM  dispersion.  A 
similar  situation  is  observed  for  the  TE  dispersion  in  (c),  as 
v>m  = 27^L_  falls  in  the  zeroth-order  gap  region.  Solid  and 
dashed  lines  correspond  to  0 = 7r/12  and  0 = 0 incidence 
angles,  respectively. 


Moreover,  if  n2B  < 0,  © is  still  valid,  with  Qa  still  given 
by  0,  but 

Qb  = ~J n2BI  - n\  sin2  6,  (11) 

where  tibi  is  the  imaginary  part  of  ns- 

For  the  TM  polarization,  the  transcendental  equations 
© and  ©,  as  well  as  the  definitions  of  Qa  and  Qb  for 
the  different  cases  are  also  valid,  provided  one  replaces  /ia 
by  €a  and  ps  by  €b,  where  ea  and  es  are  the  dielectric 
permissivities  in  layers  A and  B,  respectively. 

Figure  [2]  shows  the  calculated  dispersions  (v  = iff-)  in 
the  cases  of  TE  [Figs.  0a)  and  0c)]  and  TM  [Figs.  0b) 
anc!0d)]  polarizations,  for  different  layer  widths.  Note 
that  ve  and  vm  fall  outside  the  (n(iv))  = 0 zeroth  or- 
der bandgap.  For  9 = 0 (dashed  curves),  the  second 
photonic  band,  just  above  the  (n(w))  = 0 zeroth  or- 
der gap,  is  a pure  photonic  mode  and,  for  9 = 7r/12, 
there  is  a coupling  between  the  radiation  field  and  a plas- 
mon mode  which  leads  to  a pair  of  plasmon-polariton 
modes.  In  other  words,  while  for  normal  incidence  the 
band  just  above  the  zeroth  order  gap  is  a pure  photonic 
mode,  for  oblique  incidence  ( 9 0)  the  electromagnetic 

field  + plasmon  interaction  leads  to  a pair  of  coupled 


modes.  We  emphasize  that,  for  0 = 0,  the  electric-  and 
magnetic-field  components  along  the  z growth  direction 
are  null  and  there  is  no  longitudinal  wave  propagation. 
On  the  other  hand,  for  9 ^ 0,  in  the  case  of  the  TM 
(TE)  configuration,  the  electric  (magnetic)  field  has  a 
component  on  the  z direction,  which  results  in  the  ex- 
citation of  longitudinal  electric  (magnetic)  waves.  For 
9 = 7t/12,  Figure  0 illustrates  that  the  pair  of  plasmon- 
polariton  bands  are  asymptotic  to  the  v = ve  and  v = vm 
pure  electric/magnetic  plasmon  values.  It  is  clear  from 
Fig.  0 (a),  for  example,  that  at  small  values  of  k,  the 
lowest  9 = 7t/12  plasmon-polariton  mode  behaves  like 
an  electromagnetic  photonic  wave.  As  k increases,  the 
dispersion  bends  over  and  reaches  the  limiting  value  of 
the  pure  magnetic  plasmon.  On  the  other  hand,  the  sec- 
ond 9 = 7r/12  plasmon-polariton  branch  behaves  like  a 
magnetic  plasmon  wave,  at  small  fc,  and  like  an  electro- 
magnetic photonic  mode  as  k increases.  The  same  com- 
ments are  valid  for  the  other  panels  in  Fig.  [2]  [of  course, 
in  the  cases  of  Figs.  0b)  and0d),  the  electric  plasmon 
mode  is  the  one  involved] . By  comparing  the  results  for 
normal  and  oblique  incidence,  it  is  clear  that,  for  0 =£  0, 
resonant  plasmon-polariton  waves  are  excited  by  the  cou- 
pling of  the  electric  (or  magnetic)  plasmon  modes  with 
the  incident  eletromagnetic  field.  One  may  infer  that 


4 


the  plasmon-polariton  waves  are  excited  by  the  magnetic 
field,  in  the  TE  case,  or  driven  by  the  electric  field,  in 
the  TM  case.  This  is  consistent  with  the  fact  that  those 
longitudinal  waves,  in  the  TE  case,  are  asymptotic  to  the 
iym  pure  magnetic  plasmon  frequency  whereas,  in  the  TM 
case,  the  plasmon-polariton  bands  are  asymptotic  to  the 
ve  pure  electric  plasmon  mode. 

The  calculated  dispersions,  for  TE  modes,  as  a function 
of  the  incidence  angles  are  displayed  in  Fig.  [3]  We  no- 
tice that  the  essential  features  persist  irrespective  of  the 
incidence  angle,  although  the  Bragg  gaps  widens  consid- 
erably for  higher  incidence  angles.  Also,  Fig- 0] clearly  in- 
dicates that,  by  choosing  the  resonant  plasmon  frequency 
within  the  (n(u>))  = 0 zeroth  order  bandgap,  the  coupling 
of  light  to  plasmon  modes,  for  0^0,  essentially  disap- 
pears leading  to  a basically  pure  (electric  or  magnetic) 
plasmon  mode.  This  is  a consequence  of  the  fact  that 
the  energy  of  the  incident  electromagnetic  wave  lies  in  a 
forbidden  energy  gap  region  and,  therefore,  the  coupling 
of  the  incident  light  with  plasmons  is  expected  to  be  quite 
weak. 

Summing  up,  we  have  studied  light  propagation 
through  a ID  superlattice  composed  of  alternate  lay- 
ers of  a positive  constant  material  and  a negative  dis- 
persive material,  and  verified  the  appearance  of  cou- 
pled plasmon-polariton  modes  of  electric  and  magnetic 
natures.  The  photonic  dispersion  was  calculated  from 
Maxwell’s  equations  using  the  transfer-matrix  approach, 
with  the  LHM  modeled  by  Drude-like  dielectric  and  mag- 
netic responses.  In  the  case  of  oblique  9 =£  0 inci- 
dence, the  photonic  dispersion  indicates  that  the  mag- 
netic or  electric  field  in  the  frequency  region  around  the 


(n(u)))  = 0 zeroth  order  bandgap  leads  to  coupled  mag- 
netic or  electric  plasmon-polariton  modes  for  the  TE 
and  TM  configurations,  respectively.  Moreover,  present 
results  show  that  the  coupling  of  light  with  plasmons 
is  weakened  by  choosing  the  plasma  frequency  inside 
the  zeroth  order  gap.  As  light  propagation  is  forbid- 
den in  that  particular  gap-frequency  region,  the  coupled 
plasmon-polariton  mode  becomes  essentially  a pure  plas- 
mon mode.  This  feature  permits  one  to  select  which  type 
of  plasmon  polariton  (i.e. , electric  or  magnetic)  one  is 
willing  to  excite,  by  choosing  the  magnetic  or  electric 
plasmon  frequency  within  the  zeroth  order  gap.  Here, 
one  might  conjecture  that  the  possibility  of  excitation  of 
magnetic  plasmon  polaritons  might  provide  a new  avenue 
to  the  implementation  of  novel  techniques  and  devices 
based  on  the  interplay  of  photonics  and  magnetism  at 
the  nanoscale,  similarly  to  the  interplay  of  photonics  and 
electronics.  Therefore,  the  implementation  of  plasmonic 
chips  for  high  data  rate  processing  or  efficient  sensing  ap- 
plications might  also  be  obtained  based  on  the  excitation 
of  magnetic  plasmon  polaritons.  Finally,  we  have  studied 
an  ideal  system  in  which  losses  have  been  neglected,  and 
future  work  including  such  effects  is  certainly  forthcom- 
ing. 
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